ON MULTIPLE FREQUENCY POWER DENSITY 
MEASUREMENTS 



GIOVANNI S. ALBERTI 

Abstract. We shall give a priori conditions on the illuminations iff such that 
the solutions to the Helmholtz equation 

-div(a Vu l ) — kqu 1 = in Q, 
i 1 = ifi on (90, 

and their gradients satisfy certain non-zero and linear independence proper- 
ties inside the domain Q, provided that a finite number of frequencies k are 
chosen in a fixed range. These conditions are independent of the coefficients, 
in contrast to the illuminations classically constructed by means of complex 
geometric optics solutions. This theory finds applications in several hybrid 
problems, where unknown parameters have to be imaged from internal power 
density measurements. As an example, we discuss the microwave imaging 
by ultrasound deformation technique, for which we prove new reconstruction 
formulae. 



1. Introduction 

Let d = 2 or d = 3 be the dimension of the ambient space, il C M. d be a smooth 
bounded domain and consider the Helmholtz equation 



— div(aVit^) — kqut — in 0, 
ip on 90, 



(1) 

l k 

where a is a real and symmetric tensor satisfying the ellipticity condition and q > 
0. Let K m i a < A'max be the bounds for the possible values of k and set K a d = 
[K min , if max]- Problem ([T]) is well-posed provided that k ^ E, the set of the Dirichlet 
eigenvalues. 

We want to find suitable illuminations ip\ such that the corresponding solutions 
to |T]) and their gradients satisfy certain non-zero properties inside the domain, e.g. 

(2) l«H*)l>0, |V< 2 (x)|>0, 

and certain linear independence properties, e.g. 



(3) dct [Vu% 2 ■ ■ ■ Vw^ d+1 ] (x) > 0, dot 



«r • • • < d+i 
Vu r • • • v «r +i 



(x) > 0. 



The classical way to tackle this problem is by means of the so called complex 
geometric optics solutions [211 [10] . These are particular high oscillatory solutions of 
the Helmholtz equation and can be used to determine suitable illuminations j9[ [8] . 
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However, they can only be constructed when the parameters are sufficiently smooth. 
Furthermore, and most importantly, their construction depends on the coefficients 
a and q. Thus, this cannot be considered a completely satisfactory answer to this 
issue: in inverse problems these are usually unknown. 

The main contribution of this work is a very different approach to this problem 
by means of a multi-frequency approach. We shall give a priori, i.e. independent of 
a and q, conditions on the illuminations whose corresponding solutions satisfy the 
required properties, provided a finite number of frequencies are used in the range 
K ac i- For example, we shall show that if Q is convex then the choice {l,x 1,2:2} 
is sufficient in two dimensions, and the same is true for {1, x±, X2, X3} in three 
dimensions, provided that a is close to a constant. The main idea behind this 
method is simple: if the illuminations are suitably chosen then the zero level sets of 
functionals depending on uf, and Vu^ move when the frequency changes. Consider 
as a test case the requirement |uf(a;)| > 0. It is much easier to study the k = 
case, that is easily achieved by Maximum Principle provided that ip > 0. Then by 
using the analyticity of the map k 1— > we shall show that the required property 
transfers to any range of frequencies: it is enough to choose a finite number of 
k in [K m i n , K max \. Note that an infinite number of frequencies and analyticity 
properties have been used in similar contexts [IJ 113) , but here only a finite number 
of k is needed. Simplified versions of the main results can be found in Subsection 

ED 

The theory presented in this paper finds its applications in several hybrid imaging 
techniques. By combining measurements coming from two different modalities it is 
possible to obtain high-resolution and high-contrast images. The reader is referred 
to the recent works by Bal [8] and Kuchment [3T] for a review of the state of the 
art in hybrid techniques. Generally, a hybrid problem involves two steps. First, 
internal energies are measured inside the domain and, second, from their knowledge 
the unknown parameters have to be reconstructed. 

Many hybrid problems are governed by the Helmholtz equation ([lj or by one of 
its variants (complex coefficients, different type of boundary conditions, for which 
a multi- frequency approach can be carried out in the same way), e.g. microwave 
imaging by ultrasound deformation [25l [4] , thermo-acoustic imaging [TOl [6] , tran- 
sient elastography and magnetic resonance elastography [12] . The internal energies 
are always linear or quadratic functionals of uf. and of Vu^ and the parameters a 
and q have to be imaged. From the reconstruction procedures discussed in these 
papers, it turns out that some or all the conditions @ and ([3]) are necessary. Thus, 
being able to determine suitable illuminations independently of the unknown pa- 
rameters is fundamental, and these can be given by the multi-frequency approach 
developed here. 

As an example, we apply the theory to one of these hybrid problems, microwave 
imaging by ultrasound deformation, which was introduced in [3]. The internal 
energies have the form 

E = a\Vul\\ e = q{ul)\ 

and the electromagnetic parameters a and q have to be reconstructed. We provide 
reconstruction formulae for a/q and q, which are applicable if ([2]) and a weaker 
version of (j3]) are satisfied for a suitable set of illuminations and a finite number of 
frequencies in the microwave regime. 
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We believe that the multi-frequency approach developed in this work can be used 
in other contexts where conditions similar to ([2]) and ([3]) naturally arise [12]. In 
many situations one looks for solutions of the Helmholtz equation satisfying certain 
properties: complex geometric optics is the only tool that has been used so far [8] 
and presents the difficulties described above. 

This paper is organized as follows. In Section [2] we precisely describe the set- 
ting and state the main theoretical results. Then, we apply these to a particular 
hybrid problem and provide reconstruction formulae. In Section [3] we consider the 
mathematical aspects of the Helmholtz equation with complex k: we review well- 
posedness and regularity results and show the analyticity of the map k i— >• In 
Section |4] the multi-frequency approach is discussed and used to prove the main 
theorems. Section [5] is devoted to the study of the example and to the proof of the 
reconstruction formulae. 

Let us clarify some notation. The space of infinitely differentiable functions 
compactly supported in will be denoted by X>(fi). We use the classical notation 
for function spaces and, unless otherwise stated, we always consider real-valued 
function spaces. Whenever complex-valued functions have to be taken, we add the 
letter C, as for instance in H 1 (f2; C). Let 7 denote the trace operator from H 1 (0; C) 
to H^iQC). In a Hilbert space H, the scalar product will be denoted by ( , )h- 
We shall use the notation 6 UfV for the angle between two non-zero vectors u, v in 
a Hilbert space, i.e. 9 U>V = arccos -We write B{x,r) for the ball centered in 

x £ M. d of radius r > 0. 



2. Main Results 

Let d — 2 or 3 be the dimension of the ambient space, n 1 C n C WL d be two C 1,a 
bounded domains for some a £ (0, 1). Let a be a real, symmetric tensor satisfying 
the ellipticity condition 

(4) A|£| 2 <£-<<A|£| 2 , £ g M d , 

for some A, A > 0. We also suppose that a is of class C°- a {Tl; M. dxd ). Let q g L°°(D,) 
satisfy 

(5) /?i < q < 02 almost everywhere in fi, 

for some ^2 > 0- Let S denote the set of the Dirichlet eigenvalues of the problem 
„ f -div(a Vu£) -kqul = in fi, 

which is a countable set of positive numbers going to infinity. For any k G C \ E 
problem <jTJ) is well-posed and uf, £ C 1 (fJ;C) (see Section 3). Let ^min < K max be 
the bounds for the possible wavenumbers k and denote K a d = [K min , K max ]. More 
generally, we can assume that K a d is a non-empty continuous curve in C, containing 
the admissible values for the frequencies k and depending on the particular problem 
under consideration. 

In this section we summarize the main results of this paper. In Subsection 12. 11 
the theoretical results regarding non-zero properties of solutions of the Helmholtz 
are discussed. In Subsection 12. 21 we apply the general theory to a particular case. 
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(6) 



2.1. On the Critical Points of the Solutions of the Helmholtz Equation for 
Multiple Frequencies. Given finite subsets K C K a d \ £ and / C 7 (C 1 '" (£!)), 
the Cartesian product K x / is called a set of measurements. 

Given a set of measurements K x {tpi\, we consider the unique solution u\ € 
C 1 (il) to the problem 

— div(a Vu'j.) — kqu\ — in fl, 
u\ = (pi on dfl. 

The introduction of the following class of sets of measurements is motivated by 
several hybrid techniques [251 ITU1 SJ In this work we will not apply the full 
generality of this concept (see Definition 12.51) . 

Definition 2.1. Let p, r, s be three positive constants. A set of measurements 
K x {tfi : i = 1, . . . , d+ 1} is complete in 0' if for every x € fl' there exists k = 
k(x) £ K such that: 

(CSM 1) \u\{x)\>p, 

(CSM 2) |dct [V«| ■ ■ ■ Viif 1 ] (»| > r, 



(CSM 3) Idet 



u\ ■■■ uj +1 

k k 

Vui ••• Vu d+1 



(x)\ > s. 



k k 

The conditions (|CSM 1|) and (|CSM 2\ guarantee that meaningful power density 
measurements exist in every point. Further, (|CSM 2[) and ([CSM 3[) imply that 
there are d and d + 1 independent measurements everywhere, respectively. 

The existence of complete sets of measurements is a non trivial problem for 
general a and q. In the references cited above, some or all these conditions are 
shown by using complex geometric optics solutions |24] . In view of Proposition 
3.3 in [llj . suitable illuminations can be constructed with a fixed frequency ko, 
provided that a and q are smooth enough. However, their construction depends 
on the knowledge of a and q, that in the inverse problems we have in mind are 
unknown. 

In order to tackle this issue, in SectionUwe discuss a multiple frequency approach 
to construct complete sets of measurements. The conditions on the illuminations 
are independent of the coefficients. Here we provide a simplified version of the main 
results therein discussed. 

The next theorem discusses the construction of complete sets of measurements 
in dimension d = 2. It is a particular case of Theorem 14.81 

Theorem 2.2. Suppose d = 2, a € C 0,1 (fl) and that Q, is convex andC 2 . Then we 
can choose a finite K C K a d \ £ such that 

K x {l,x-y,x 2 } 

is a complete set of measurements in fi. 

The next theorem deals with the construction of complete sets of measurements 
in dimension d = 3. It is a particular case of Theorem 14. 121 

Theorem 2.3. Suppose d = 3. If a is constant then we can choose a finite K C 
K a d \ S such that 

K x {l,xi,x 2 ,a;3} 
is a complete set of measurements in Q. 
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The proofs of these statements can be found in Section |4j which also contains 
more general versions of the previous results. In particular, we shall see that in two 
dimensions the illuminations l,x\ and Xi represent a simple selection from a much 
wider class and the domain does not have to be convex. In the three-dimensional 
case, a is in fact required to be close to a constant matrix. Moreover, we describe 
how to choose the finite set of frequencies K. 

2.2. Applications to Microwave Imaging by Ultrasound Deformation. We 

now discuss an example to justify the introduction of complete sets. In fact, for 
our purposes, a wider class of sets of measurements is sufficient. 

We consider the hybrid problem arising from the combination of microwaves and 
ultrasounds which was introduced in [4] . Full details and the proofs of the results 
are given in Section[5j In addition to the previous assumptions, we suppose that o is 
scalar- valued. In microwave imaging, a is the inverse of the magnetic permeability, 
q is the electric permittivity and K ac i = [if m i n , if max], with K m i n > 0, are the 
admissible frequencies in the microwave regime. 

Given a set of measurements K x {ifi} we consider power density measurements 
of the form 

(7) ^ = qului, E^aV^-Vni, 

where u\ is given by ([6]). For simplicity, we denote = and similarly for E. 

These internal energies have to be considered as known functions in fi'. In contrast 
to [3] , we do not assume the cross- frequency data to be available: thus the available 
data is smaller, and this complicates the reconstruction. 

Problem 2.4. Choose a suitable set of measurements K x {(fi} and find a and q 
in fi' from the knowledge of e l i and EV in fi'. 

Problem l2.4l is solved via two reconstruction formulae for a/q and q, respectively, 
which we shall now describe. Their applicability is guaranteed if K x {fi} is a proper 
set of measurements in fi'. 

Definition 2.5. Let p and r be two positive constants. A set of measurements 
K x {ipi : i = 1, 2, 3} is proper in fi' if for every a; € fi' there exists k = k(x) £ K 
such that: 

(PSM 1) \u\(x)\>p, 



(PSM 2) |Vu|(a;)||Vu|(a;)| sin 6» Vu 2 v „ a (x) 



> r. 



The collection of all proper sets of measurements in fi' with constants p and r will 
be denoted by V(fl';p, r). 

Clearly, the conditions characterizing a proper set are weaker than the conditions 
of a complete set: the main difference relies in the requirement of only two inde- 
pendent measurements in any dimension. Thus, the construction of proper sets can 
be easily achieved by using the Theorems of the previous subsection (in dimension 
three the illumination x% is not needed). 

The next statement gives an exact formula for a/q. The formula was first derived 
in [4] for the two-dimensional case, and we have extended it to any dimension. We 
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use the following notation for a matrix-valued function M of size N 

N N 

tr(M) = ]T M,,, \VM\ 2 2 = ]T | VM y | 2 . 
Theorem 2.6. Take p,r > and K x {(pi} G V(^l';p,r). Suppose that 

(8) |if (z)!,!^)! <&, iGfy, 

for some b > 0. Tafce x € O' and k as in Definition \2.5l Then there exists 
C = C(X,/3i,b) > such that 

(9) tr( e ,)tr(£,)-tr(e^) ( ^ cA , 

tr ( e fc) 

and a/q is given in terms of the data by 



,2 a tr(e- k )tr(E- k )- tr(e h E n ) 

q tr(e k y 



(10) | V (e 8 /tr(e^)|^=2 ^^^^ m x. 



We now give an exact formula for q. If compared to the one described in [3], this 
one is valid in any dimension and does not require the set of measurements to be 
complete. It only requires that 

tr(e) := eji? > c> in 

i,h 

which holds provided that the set of measurements satisfies (jPSM II) . 

Theorem 2.7. Let K x {tpi} be a proper set of measurements. Suppose q £ H 1 (f2) . 
T/ien log g is the unique solution to the problem 

-div(Gtr(e) Vu) = -div (G V (tr(e))) + 2J2 k d (E% - ke\ l ) in Q', 
u = logman' on dQ,'. 

The previous result allows to reconstruct q from the knowledge of G = a/q. Also, 
<? is supposed to be known on the boundary of f2' , which is a reasonable assumption 
since we are mainly interested in detecting inclusions inside the domain. At this 
point, an obvious formula for a is a = Gq. 

The two exact formulae provided by the previous Theorems use the derivatives 
of the data. As a consequence, in presence of noise, they may not give good qual- 
ity images. Thus, an optimal control approach would be an efficient way to find 
better approximations of the coefficients, starting from the good guess given by the 
formulae described here. Details can be found in pQ. 



3. The Helmholtz Equation 

In this section the Helmholtz equation ([1) is discussed, where a £ C°' a (£l;M. dxd ) 
and q £ L°°(f2) are as in Section [2] and satisfy ((4]) and <|5j) . First, we show that 
the problem is well-posed in the general situation of complex k. Then, regularity 
results are discussed. In order to do this, we use the classical results on elliptic 
problems. 
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3.1. Well-Posedness and Regularity. First of all, we study existence, unique- 
ness and stability with homogeneous Dirichlet boundary conditions. The space 
i? _1 (f2; C) denotes the continuous antidual of Hq(Q.; C). 

Proposition 3.1. There exists £ = {Xi : i G N} C R + with A; — > +oo such that 
for fc G C\S and f G ff _1 (fi;C) the equation 

(11) — div(aVu) — k qu = f 

has a unique solution u G Hq(Q.\C) satisfying \\u\\ Hl < C\\f\\ H -i for some C > 
independent of f . Moreover, for fixed f G C), the map 

fc G C\ £ i — >u G H^{Q;C) 

is analytic. 

Proof. We give a sketch of the proof. Take L = -div(oV-): H^(Q:C) ~> 
and Mq : L 2 (fi; C) -> L 2 (fi; C) defined by / i-> qf. By Lax-Milgram Theorem and 
(|4]) the operator L is invertible and we can consider S := L~ 1 AI q : Hq(&1;C) — > 
i?o(fi;C), which by Kondrachov Compactness Theorem and ([5]) is compact, self- 
adjoint and positive. Hence S has a countable set of eigenvalues {rji > : i G N}, 
with r)i -> 0. Define S = {A, = 1/r?,- : i G N} and take HC\E. Since dTTJ) is 
equivalent to (I — kS)u = L^ 1 f . the first part follows. Finally, the analyticity 
of k n- u is a consequence of the so called Analytic Fredholm Theorem (see |23[ 
Theorem 8.92]). □ 

As a consequence, we immediately get the following result regarding the Dirichlet 
boundary value problem for the Helmholtz equation. 

Corollary 3.2. Take k G C \ S. Then there exists C = C(fi, a, q, k) > such that 
for every f G _ff _1 (fi; C) and tp G -ff 1//2 (f2; C) the problem 

J — div(a Vu) — k q u = f in£l, 
\ u = Lp on d£l, 

has a unique solution u G i/ 1 (r2;C) satisfying 

(13) IMIffi(Q;C) - C (lMlffi/2(fi ;C ) + l/llff-i(f2;C)) • 

Moreover, for fixed f G C) and ^ G H 1/2 (fl; C), i/ie map 

fc g c\s i — >ue ^(njC) 

is analytic. 

Standard elliptic regularity theory allows us to study the regularity of the solu- 
tion u G iT^fijC) to (fl2|) . 



Proposition 3.3. Tafce fc G C \ S, / G L°°(fi;C), u G C 1,a (fi;C) and rante 
V = l( v )- Let u G i? 1 (f2; C) oe i/ie unique solution to ilty) . Then u G C 1 (fi; C) and 
£/iere exists C = C(f2, a, g, fc) > smc/i i/ia£ 

(14) ll«ll C i(n iC )<c(lMI C i.«(n :C ) + ll/H 
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Proof. We decompose u, /, tp and ft into real and imaginary parts by writing u = 
ur + i ui, f = fn + i fi, v = vr + ivi and ft = Ur + i fc/. By testing the first 
equation of (|12[) against any w G Hq(Q; R) and by taking real and imaginary parts 
we obtain 

\ -div(aVufl) - k R qu R = -kiquj + / R , 7(ur) = 
|-div(aVw/) - k R qui = kj qu R + /j, 7(7/7) = 95/. 

In the rest of the proof we will consider constants c = c(fi,a, g, ft). From usual 
elliptic regularity theory (e.g. (8.38) in [TS] ) there holds 

IMloo ^ c (JMI2 + \\ v \\c^(n;C) + ll/lloo) 
By arguing in the same way with ui we infer that 

IMloo <c(H 2 + IMI C i, Q( n ;C ) + 11/11 

Therefore, by Corollary 8.35 in [H5] we get that ur, uj G C^fi) and that 

\\u R \\ c i m < c (\\u\\ 2 + \\v\\ cl . a(n . c) + ll/IU). 
A similar inequality holds for uj and so 

IMI C i(fi ; C) ^ C (NI2 + IMIc±.«(Q;C) + ll/lloo) ■ 

Finally, in view of {13]) we obtain (jT4j) . □ 

3.2. Analyticity Properties. First, we need the following lemma concerning the 
composition of analytic functions. 

Lemma 3.4 (|27j). Let D C C be open, X,Y be Banach spaces, /j : D — > X be 
analytic maps for i = 1, . . . , b and g: X b — > Y be multilinear and bounded. Then 

.go (/!,...,/„): D^Y 

is analytic. 

In this subsection we study the dependence of the solutions u\ on k. We come 
back to the original problem 



div(a Viij) — kqu^ = in fi, 



= (p on <9fi, 



for ft G C \ E and fixed </? G 7 (C^fi)). By Corollary 0H there exists a unique 
solution G i/ 1 (fi;C), that in view of Proposition 13.31 is in C 1 (fi;C). We have 
already seen that uf depends analytically on the wavenumber ft with respect to 
the norm of if 1 (fi;C) . We now want to show this with respect to the C : (fi;C) 
norm. The proof follows the argument used by Calderon to prove that u\ depends 
analytically on a if ft = |20| . 

Proposition 3.5. Take tp G 7(C 1 ' Q (fi)). Then the map 

T:C\E — ^(fijC), k 1 — >ul 
is analytic, where is the unique solution to {Ip. 
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Proof. For simplicity, we do not write the dependence on (p. Fix ^ £ C \ E: we 
shall prove that T is analytic in feo- 

For h G C define the operator C h : C^HjC) L°°(fi;C), w ^ hqw with 
norm < /?2 H- In view of Proposition 13.31 there exists a bounded operator 

B: L°°(£l; C) — > C 1 (H; C), where -B/ is the unique solution to 

— div(a Vu) — k qu = f in f2, 
u = on 9fJ. 

As a consequence, BC\: C 1 (fi;C) — > C 1 (fi;C) is a bounded operator with norm 
\\BC h \\ < fo\\B\\\h\. Define r = min {l/(/3 2 ||B||), d(k , £)}, where d(*o,£) = 
inf {c?(fco, a) : a G E} denotes the distance between fco and E. Take k G B(ko,r) C 
C \ E. The increment h = k — ko satisfies < r < ^ ^gri , whence H-BC^H < 1. 
Denoting v = — uu a , by definition of T we have 

-div (a Vw) — koqv — hqv = hquk = C ft it fco 

which, after applying B to both sides, becomes (7— BCh)v = BChUk , where I 
denotes the identity operator of C 1 (fi;C). As a result, since H-BC^H < 1, we obtain 

oo 

v = (I - BC h y l BC h u ko = ( BC hT u k«- 

n=l 

It follows that 

oo 

(15) r(fc) = ^(BC^ D )"r(fco), keB(k ,r). 

71=0 

Namely, the map T is analytic. □ 

For the sake of completeness, we now show that the dependence on k of the 
internal energies ([7]) is analytic. Consider two different boundary data ip\,if2 G 
_ff 1 / 2 (ri;C) and the corresponding solutions u\ := u^' to ([I}. For i,j = 1,2 we 
define the internal energies as 

e 13 : C \ E — > L x (0; C), E ij : C \ E — > L 1 (f2; C), 

k i — > q u\ ttl, fc i — > a Vitj. • VuJ. 

Theorem 3.6. Tafce ipi,(f2 G iJ 1 / 2 (J7;C). TTien £/ie internal energies e 4 -? and E % i 
are analytic functions. 

Proof. Note that the maps 

3i : H 1 ^; C) 2 — -> L x (0; C), g 2 : ff 1 ^; C) 2 — > L 1 ^; C) 

(u, w) i — )■ g tt v (it, v) i — > a Vu • Vit 



are bilinear and bounded. Therefore the result follows by Corollary [321 and Lemma 

□ 



As a consequence of the unique continuation property for holomorphic functions 
we obtain the following 

Corollary 3.7. Take <p±,(f2 G H 1/2 (VL\ C). Let ki, k G C \ E such that ki k and 
ki =/= k for every I £ N. If e lJ [resp. E lJ j is known in ki for every I G N then e lJ 
[resp. E lJ ] is known everywhere in C\ E. 
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Remark 3.8. This result must be seen in view of Problem 12.41 The knowledge of 
E lJ for infinitely many frequency in a fixed range determines E lJ for k = 0, where 
the reconstruction process for a has been studied thoroughly [5J [51 [HI El [221 HE] • 
However, this has to be regarded simply as an interesting theoretical result: analytic 
continuation is a very ill-posed process. 

4. On the Critical Points of the Solutions of the Helmholtz 
Equation for Multiple Frequencies 

In this section we discuss the multi-frequency approach to the issue of the ex- 
istence of complete and proper sets of measurements. The theorems stated in 
Subsection 12.11 will be a consequence of the results presented here. 

4.1. Preliminaries. Let a £ C°- Q (fi; R dxd ) and q £ L°°(ft) be as in Section [2] 
and satisfying ([4]) and ([5]). Recall that K a d is the set for the possible wavenumber, 
namely we are allowed to choose values k £ i^ a[ j\S. As in Proposition ^. 51 we denote 
the unique solution to the boundary value problem (p} by u\ = T v (fc) £ C 1 (fi; C) 

„ f -div(aVw^) - kqul = in 0, 

As a consequence of the analyticity of k i— > proved in Proposition 13.51 we 
obtain the following 

Lemma 4.1. Take Q' C fi, 6 £ N* and (: C^HjC) C (fFjC) multilinear and 
bounded. Let <pi, . . . , (fib £7 (C 1,Q (f2)) be such that 

(16) C «P) (*)^o, xeW, 

for some P6C\S. Tafce fc;, fc £ C \ S with ki — > fc and fc; 7^ fc. TTien there exists 
a finite LCff such that 

(17) EK«V •• '00*01 >0, xeT?. 
leL 

In particular, we can choose a finite K C K a d \ S such that 

(18) £|C«\...,< i ')(*)|>0, xeW. 

fee A' 

Remark 4.2. For the sake of completeness, we note that the statement holds true 
for £ analytic, but our current applications do not need this extension. We shall 
use this Lemma only with k x = 0. However, this more general version allow other 
possible applications. In particular, the case \k x \ — > 00 could be studied. This could 
be considered in the context of high-frequency approximations of the Helmholtz 
equation. 

Remark 4.3. If in addition to the assumptions of this lemma we assume that the 
set {x € f2' : ( (u^ 1 , . . ■ , u^ b ) (x) = 0} is finite, then there exists I £ N such that 

|C (up ,...,<) (*)l + |C(< >.•• , u£) (x) I > 0, X £ W. 

Namely, if the zeros of ( (u^ 1 , . . . , u^ b ) are isolated points then only two frequencies 
are necessary to obtain a non-zero quantity everywhere. We leave the proof of this 
extension to the reader. 
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Proof. Take x £ f2'. The map ev x : C ffi'; C) — s- C, u n- is linear and bounded. 
Hence, using Proposition 13.51 and Lemma 13.41 the map 

z x :=ev x o(o(T^\...,T^) : C\E^C 

is analytic. By (jTSJ) we have z x (k x ) ^ 0, and so the set {k € C \ £ : ^(fc) = 0} has 
no accumulation points in C \ E. As a consequence, there exists S N such that 
z x{ki x ) 7^ 0. Since ( o (T Vl , . . . ,T Vb ) (fc; x ) is continuous, we can find r x > such 
that 

(19) z y (kj^0, yeB(x,r x )r\W. 

Since ft' = U^gTF ^ there exist a;i, . . . , xn £ ft' satisfying 

N 

(20) TF=\J(B{x l7 r Xt )nTF) . 

1=1 

Defining L = : i = 1 : iV} . by ([15]) and ([2"0f we obtain (JT7]). □ 

Let us now recall the definition of complete sets of measurements. 

Definition 12.11 Let p, r, s be three positive constants. A set of measurements 
K x {ifi : i = 1, . . . , d + 1} is complete in 0' if for every x £ il' there exists k = 
k(x) € K such that: 

(CSM1) \ul(x)\>p, 

(CSM 2) |det [Vuf • • • Vw^ +1 ] (x)\ > r, 

(CSM 3) |det 

In order to satisfy these conditions, the main idea is to trace back to the case 
k = 0, where things are simpler. For instance, consider condition (jCSM lj) . With 
k = 0, by the Strong Maximum Principle this is trivially satisfied provided that 
the boundary condition is strictly positive or negative. It remains to show that 
this good property transfers to any range of frequencies, and this is achieved with 
Lemma l4.ll In summary, the strategy to study these conditions can be outlined in 
the following three steps: 

(1) choose a suitable £ as above such that the condition we want to prove is 
equivalent to (fT5| : 

(2) prove the assumption (|16p for k x = 0, which is in general easier than (| 18|) 
since k = 0; 

(3) apply Lemma |4~T1 to deduce the result. 

We will deal with the issue of the construction of complete sets of measurements 
in two different situations, depending on the dimension. As we shall see, condition 
(|CSM 1|) will be a consequence of the Maximum Principle, that holds in any di- 
mension and for any a. However, the study of conditions ([CSM 2j) and ()CSM 3[) in 
the case k = depends on the dimension. 



4 



.d+i 



d+l 



(x)\ > S. 
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4.2. Complete Sets: Two-Dimensional Case. Throughout this subsection we 
assume d = 2. First, we consider conditions (|CSM lj) . (|CSM 2j) and (|CSM 3|) for 

= 0. As far as (jCSM 1[) is concerned, the Maximum Principle will be the main 
tool, and no further investigation is required. 

Let us now focus on (jCSM 2|) for k = 0, which reads 

det [Vu 2 , Vug] > r 

for some r > 0. This problem has been studied by Bauman et al. [M] in the context 
of univalent mappings. We introduce the following class of boundary conditions. 

Definition 4.4. Let V = (^2,^3) G C 1,Q (0;R 2 ) be a C 1 diffeomorphism of a 
neighborhood of fl into R 2 such that Jip > in 17. Denote the restriction of i/j to 
the boundary dfl by (p — {(f2, ^3) := il>\dn- We say that ^2, ¥>3 are BMN boundary 
conditions if ip: dfl — > R 2 maps 90 onto a convex closed curve. In this case we 
write (^2,^3) G BMN(Q). 

Remark 4.5. Note that if O is a convex domain then tp*2 = %i, tp 3 = £2 are BMN 
boundary conditions. 

Their main result shows that (|CSM 2|) is satisfied for k = 0, provided that the 
boundary conditions are BMN. 

Proposition 4.6 (Theorems 2.4-2.5, Q3). Suppose d = 2. If {ip 2 , ip 3 ) G BMN{Vt) 
then 

det [Vug Vug] > in fi. 
// m addition Q is a C 2 domain and a e C ' 1 ^; R dxd ) i/ien 

m fi. 



det [Vu 2 . 



Vug] > 



Let us now focus on (jCSM 3[) for k — 0. As far as the author is aware, there is 
no result concerning the zeros of 



det 



Vu 2 



Vug 



Therefore we need the following statement, whose proof can be found at the end of 
this subsection. 

Lemma 4.7. Suppose d = 2. Take ipi,tp2,<f3 <G 7 (C l a (Q)) such that ipi has a 
fixed sign and (p 2 , <p 3 ) , G BMN(Ct). Then 



dCt [VuJ Vu 2 VugJ ^ U m "- 

// m addition is a C 2 domain and a G C ' 1 ^; R dxd ) i/ien 



det 



VuJ 



Vu 2 



Vug 



O. 



We are now in a position to state the main results of this subsection, concerning 
the construction of complete sets of measurements in dimension d = 2. 

Theorem 4.8. Suppose d = 2. Take O' <s fl and ^1,^2,^3 G 7 (C 1,a (0)) such 
that tpi has a fixed sign and ((fz, <p 3 ) , G BMN(Q). We can choose a finite 

K C K a d \ £ such that 

K x {tp u (p 2 , ip 3 } 
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is a complete set of measurements in fl' . 

If in addition f2 is a C 2 domain and a G C 0,1 (S1; ]R dxd ) then we can choose a 
finite K' C K ac i \ £ such that 

K' x {</>!, <p 2i ip 3 } 
is a complete set of measurements in Q. 

Remark 4.9. Note that the hypotheses on the boundary conditions are satisfied if 
ipi = 1 and ((p2,<fi3) G BMN(il). Therefore, if f2 is convex, by Remark 14.51 an easy 
choice for the boundary conditions is tpi = 1, tf2 = x\ and ip 3 = x%- 

Proof. By Lemma 14. 7\ Proposition 14.61 and Strong Maximum Principle we obtain 
that 



4, det [Vmq 



Vug] , det 



Vug 



Vug 



^0 



o. 



respectively. Thus we can apply Lemma |4 . 1 1 with f2' <g SI, b = 6, fc x 

¥>5 = ^2, <^6 = ^3 an d C( u i w i w 7 z i v) = wdet [Vi> Viw] det 



obtain the existence of a finite set AT C K a d \ £ such that 



V/ 



= 0, <^4 =_<Pl 

z y 
Vz Vy 



to 



p(x) := det [Vi 



Vuf] det 



Vu£ 



Vu 2 fc 



(x) > 0, 



x G ST 



In the sequel, we shall denote several positive constants independent of x G SI' by 

??(x) > c for every x G 0'. As a result, for 

(x) > c. 



c. Since G C 1 (SI) (Proposition 
any x G SI' there exists fc G K such that 

4 det [Vu 2 Vuf] det 

As a consequence we have 



\ul(x) \ > c, | det [V? 



> c, 


det 







v«2 



(x) 



> C 



whence the first part of the theorem follows. 

If in addition we suppose that SI is a C 2 domain and that a G C 0,1 (S1; M. dxd ), we 
can use the second part of Proposition 14.61 and the second part of Lemma 14.71 to 
infer that 



det 



Vu 2 



Vug 



, det [Vuf t 



Vul] ^ in SI, 



respectively. Therefore, arguing as before, we obtain a complete set in O. 



□ 



From the proof of Theorem 14.81 it is clear that the assumptions ipx ^ and 
(ip 2 , <p 3 ) G BMN(ft) allow us to deduce (|CSM lj) and (|CSM 2|) of Definition [23J re- 
spectively. One may wonder if the above Theorem still holds without one of these 
two hypotheses. The answer is no, as the following examples show. Therefore, 
if the boundary condition ifi are not chosen properly, in general one cannot ex- 
pect to obtain (jCSM lj) and (jCSM 2j) by doing many measurements with different 
wavenumbers but with fixed illuminations. First, we provide a counterexample for 
condition (jCSM lj) . 
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Example 4.10. Suppose d = 2, Q = B(0, 1) and a = q = 1. We choose the 
boundary condition 991 (xi, X2) = xi, that clearly does not have a constant sign. In 
this case the Helmholtz equation ([1} can be written in polar coordinates (r, 9) and 
reads 



2 + fcu = 0. 
It is straightforward to see that 



1/ ^ Ji(Vkr) 

satisfies the above equation, where J\ is the Bessel function of the first kind of 
order 1 and k > is not an eigenvalue of the problem. Thus, {uj,} represents a 
family of solution to the Helmholtz equation with fixed boundary condition and 
varying wavenumber. However, condition (jCSM lj) cannot hold since u\(0,X2) = 
for every k and for every X2 € (—1,1). 

Next, we study condition (|CSM 2ft . Since it expresses the linear independence of 
the gradient of the solutions inside the domain, we shall see that it is not possible 
to require ipi, ip% to be just linearly independent. 

Example 4.11. We consider the situation of Example 14.101 Suppose £1 = B(0, 1) 
and a — q = 1. We choose the boundary conditions (^2(^1,2:2) = %i and (p$ = 
1. Clearly, (^2, ¥3) ^ BMN(Sl) since ip^ is a constant, but tpz, cps are linearly 
independent. It is straightforward to see that the corresponding solutions to ([1]) 
are 

Ji(Vfc) Jo(vk) 
where J„ is the Bessel function of the first kind of order n and k > is not an 
eigenvalue of the problem. 

Take a matrix-valued function A: fl — > GL(2), where GL(2) denotes the set of 
2x2 invertible matrices. By viewing A(x) as a change of coordinates in T^O, the 
tangent space in x to f2, we get 

det [AVulA- 1 AVulA- 1 } = det [A [Vu 2 k Vu\] A- 1 ) = det [Vu 2 k Vu\] . 

Therefore, as far as det [Vu^ ^ U V\ ' IS concerned, we can express the gradient in 
any system of coordinates. 

In this case, writing Vu^, with respect to eg and e r we have Vu|. = \^gg- e e + 
^ UlL e r . Hence 



dr 



det [Vul S7u 3 k ] = det 



'idul idul 




_i 


r d0 r 06 


= det 


r 



Ji(Vkr) ■ , 



Thus, we have det [Vu| Vm^] (xi,0) = for every k and for every x\ € (—1, 1) 
and so (|CSM 2[ ) cannot be satisfied by using many measurements with these fixed 
illuminations and varying wavenumbers. 

We end this subsection with the proof of Lemma 14.71 

Proof. (Lemma \4- 7| ). We shall prove only the first part of the lemma. The second 
part can be proved following the same argument. 
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For simplicity of notation we write Ui := Uq and suppose ipi > 0. By Proposition 
14.61 we have det [V«2 VU3] > in 17, and so {Vit2, VM3} form a basis of M 2 in 
every point of il. Therefore there exist A,/i: il —> R such that 

(21) Viti = AVw 2 + ^Vm 3 in ft, 

whence 



det 



Ml u 2 u 3 

Viti Vlt2 Vl*3 



(mi — Xu2 — /iU3) det [Vit2 VU3] in fi. 



By contradiction, suppose that there exists xq € fl such that 
(22) ui{x ) = \(x )u 2 (x ) + fj,(x )u 3 (x ). 

Denote Ao = X(xq), = /i(xo) and z = A0M2 + A*o M 3- Since ipi > 0, by the Strong 
Maximum Principle we get ui > in fl. Let h = zju\ and ipi ~ Ui/ui for % = 2, 3. 
For every w € 2?(S1) and i = 2,3 a straightforward calculation shows that (see 
Lemma I5.ip 



u 1 aVV'i • Vi> dx = 0. 

Therefore ipi is the unique solution to the problem 

-&i~v(u\ aWipi) = in Q, 



ipi = ifi/ifi on <9fi, 
and by using Proposition 14.61 we infer that 
(23) det [VV>2 V^] > in 17, 

since (^7,^7) € BMN(Q). By (EU) and (H2J) we have Vui(aJ ) = Vz(i ) and 

^i(a;o) = ^(^o) > 0, whence V(logMi)(a;o) = V(logz)(a;o) and so V(log h)(xo) = 0. 
As a consequence, = \7h(x ) = A Vi/'2(^o) + Mo^^^o), which contradicts (|2"3")) 
since (A ,A*o) 7^ (0,0). □ 

4.3. Complete Sets: Three-Dimensional Case. Throughout this subsection 
we assume d = 3. First, we need to consider (jCSM 2|) and (jCSM 3)) for A; = 0. 
In contrast to the case d = 2, we cannot use Proposition 14. 6[ since the equivalent 
statement is known not to hold in three dimensions |15) . 

Thus, we assume that the parameter a is a small perturbation of a constant, 
symmetric and positive definite matrix ao, i.e. 

a = a + s, \\s\\ c0 , a(u . Rdxd) < 5, 

with 6 small enough. The study of this approximation is common in the literature 
[231 HU1 [31 [5]. Note that we do not make any assumptions on g. 

In order to understand why this approximation is useful, let us consider the 
constant case, i.e. a = ciq. Choose ipi = 1, (p% = x\, (p$ = x% and tp^ = X3. Hence 
1, X\, X2 and X3 are the solutions to the problems 

— div(ao Vti) = in 0, 
u = ifi on dfl, 

for i = 1, ... ,4, respectively. Therefore conditions (jCSM ip . (|CSM 2[) and (jCSM 3j) 
are trivially satisfied in the case k = 0. Thanks to the multi-frequency approach 
we can extend this property to any range of frequencies, and a continuity argu- 
ment allows small variations around a constant value. These two steps are carried 
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out in the following theorem, which concerns the construction of complete sets of 
measurements in dimension d = 3. 

Theorem 4.12. Suppose d = 3. Let a$ be a constant, symmetric and positive 
definite matrix. There exists S > such that for any s G C°' a (Q' 1 M. dxd ) with 
\\ s \\c°- a (n-s. dxd ) — ^ we can c h° ose a finite K C K ad \ T, such that 

K x {l,xi,a;2,X3} 
is a complete set of measurements in Q for a = ao + s. 

Proof. Fix < e < A/3, where A is the smallest eigenvalue of ao. In the following 
we will consider positive constants c depending on ao, f2 and e. Take < 5 < e 
and s G C 0,a (£l;M. dxd ) with ||s|| c o, a^dxj) < S. For a s = ao + s consider u l (s) the 
solution to the problem 

J — div(a s Vuq(s)) = in f2, 
\ Uq(s) = on 9fi, 

where = 1 and ^ = 2Jj_i for i = 2,3,4. Thus, 1 1 C- 3 ) 1 1 2 — c ^ or some c > CI- 
As a result, from usual elliptic regularity theory (e.g. (8.38) in [19]) we obtain 
ll M o( s )|| 00 — c - By Theorem 8.33 in [TH] this implies that 

( 24 ) ll V «o(s)|| c o, Q(n;Rd) < C . 

The same argument applied to u' l (s) — Uq(0) as a solution to 

-div (a V(w (s) - w ( ))) = div(sVw (s)). 

gives ||mq(s) — u (0)|| < cS for some c > 0. Hence, by (f2"4"|) and Theorem 8.33 in 
[T9] applied to Uq(s) — u (0) there holds 

\\Vul(s) - V Uq (0)L < c - 4(0)1^ + 11^0^)11^^)) < cS, 

whence 

||«o(«) -«o(°)|| c i^ " C(5 ' 
As a consequence, since for s = we have 

wj(0) = det [Vug(0) • • • Vu^(0)] = dot 

there exists S < s such that 



^ in 0, 

for all s G C 0,O! (il; ]R dxc ') with 1 1 s 1 1 c0 , <=. (?7-]R<i x <i) < ^- Finally, arguing as in the proof 
of Theorem 14. 8[ the result follows by using Lemma 14.11 □ 

4.4. Proper Sets of Measurements. This subsection is devoted to the study of 
proper sets of measurements, which were introduced in Subsection 12.21 Since the 
conditions characterizing proper sets are weaker than those of complete sets, we are 
not going into the details: the reader is referred to the previous parts of this section. 
We first recall the definition of proper sets of measurements for a G C 0,a (f2; R dxd ) 
and q G L°°(0) satisfying (@} and ©. 



«S(o) 
v«5(o) 



u$(0) 
Vug(0) 



1 in Q. 



ul(s), det [Vul(s) ■■■ Vu^(s)],det 



u o( s ) 



u o( s ) 
Vug(«) 
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Definition 12.51 Let p and r be two positive constants. A set of measurements 
K x {ipi : i = 1, 2, 3} is proper in Cl' if for every x £ Cl' there exists fc = fc(x) € A' 
such that: 

(PSM1) |«je(*)| >P, 



(PSM 2) |Vu|(a;)||Vit|(a;)| sin 6» Vu 2 Vm3 (x) 



> r. 



The collection of all proper sets of measurements in Cl' with constants p and r will 
be denoted by V(Cl';p, r). 

The multi-frequency approach discussed in the last section is easily applicable 
for the construction of proper sets of measurements. In general, the choice for the 
illuminations l,xi and xi gives a proper set in any dimension (in dimension three, 
a is required to be close to a constant). 

Note that, in general, condition (|PSM 2[) is weaker than (ICSM 2|) . and they are 
equivalent when d = 2 since 

(25) |det([to z])\ = \w\ \z\ \sva.6 w>z \ , w,zeR 2 . 

Thus, when d = 2 one can directly use an analogue version of Theorem l4.8l without 
the assumption (^, ^) e BMN(Cl), which was needed to satisfy (|CSM 3|) . 

Theorem 4.13. Suppose d = 2. Tafce O' <s O and y>i, (p2,<P3 & 7(C 1,Q (n)) suc/i 
</iai <pi /ias a fixed sign and (^2,^3) £ BAIN (CI). We can choose a finite K C 
A a( j \ E smc/i f/iaf 

A X {tp u <p 2 , (f 3 } 

is a proper set of measurements in CI' . 

If in addition Cl is a C 2 domain and a G C°^(Cl; R dxd ) then we can choose a 
finite A' C A a( j \ S such that 

K' x {tp 1} tp 2 , ^3} 
zs a proper set of measurements in Cl. 

In dimension three, it is straightforward to check that the analogue version of 
Theorem 14. 121 is applicable, without the illumination z. 

Theorem 4.14. Suppose d = 3. Let ao be a constant, symmetric and positive 
definite matrix. There exists 5 > such that for any s £ C°' a (Cl;R dxd ) with 
ll s llc a (n-R t!X£i ) — ^ we can choose a finite A C K ad \ E such that 

K x {1, £1,0:2} 

is a proper set of measurements in Cl for a = ao + s. 

4.5. Numerical Experiments on the Number of Needed Frequencies. The 

theory developed so far gives conditions on the illuminations to construct complete 
and proper sets of measurements, provided a finite number of frequencies are chosen 
in a fixed range. However, we do not provide an estimate on the number of the 
required frequencies, i.e. the cardinality ffK of K. In order to document this 
point, we have performed a numerical test in 6561 different cases to see what we 
can expect. For simplicity, we have decided to consider the two-dimensional case 
and the concept of proper sets of measurements. 
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FIGURE 1. A particular combination of the coefficients. 




Take = 23(0, 1). We now describe the coefficients we have used. Set c = 0.35, 
r = 0.2 and P\ = (— c, — c), P 2 = (— c, c), P3 = (c, — c) and P 4 = (c, c). Let Xi be 
the characteristic function of B(Pi,r) (or, more precisely, a smooth approximation 
of it). Then we set 

4 4 

i=l i=l 

where e^,/?, € {0,f,2} (see Figure [T]). This construction gives 3 s = 656f different 
combinations. 

In view of Theorem 14. 131 we choose the illuminations ipi = 1, ip% = x\ and 
ips = X2 and the frequencies between the first and the second eigenvalue, in order 
to simplify the numerical computation. For any combination of the coefficients, we 
compute the number of needed frequencies to obtain a proper set of measurements 
in SI. The results are summarized in Table[TJ These figures suggest that in practical 
applications the number of needed frequency could be quite small. 



Table 1. Number of combinations of coefficients per number of 
needed frequencies to obtain a proper set of measurements in Q. 



Needed frequencies (#K) 


2 


3 


> 4 


Combinations of coefficients 


1609 


4952 






5. Applications to Microwave Imaging by Ultrasound Deformation 

This section is devoted to the discussion of the hybrid problem introduced in 
Subsection 12.21 and to the proofs of the results therein presented. 

5.1. Formulation of the Problem. The object under examination is a smooth 
bounded domain f2 C R d , for d — 2 or d — 3. In the microwave regime, the electric 
field uf, in Q is assumed to satisfy the Helmholtz equation 

, . f -div(a Vu£) - kquf, = in 0, 

^ ' \ u£ = (p on dQ, 
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where a <E C°' a (£l) is the inverse of the magnetic permeability, q € L°°(n) is the 
electric permittivity and \fk is the wavenumber (with an abuse of language we shall 
sometimes refer to k as the frequency or the wavenumber). We assume © and ([5]). 
The use of the Helmholtz equation is a very common scalar approximation of the 
Maxwell equations [TBI 121)1 [U [5] . As we have seen in Section [31 problem ([2"B|) is 
well-posed provided that k is not a resonant frequency. 

Practitioners can choose a frequency k > in a fixed range (not a resonant fre- 
quency), a real illumination ip on the boundary and measure the generated current 

on the boundary a-g^-- As described in [3], these measurements are combined with 
localized ultrasonic waves focusing in small regions oj inside fi. We assume that 
the electromagnetic parameters are affected linearly with respect to the amplitude 
of the ultrasonic perturbation, that is supposed to be small. Moreover, this mod- 
ification is localized only in the region lu. Under these assumptions, denoting the 
modified coefficients by a and q we have 

a = a (1 + c a axu>) , 
q = q (l + C q aXu) , 



where a is the amplitude of the ultrasonic signal and c a and c q are known functions. 

■ 



The corresponding electric field ut is the solution to 



— div(a V«k) — kquf, = in f2, 
u£ = ip on dfl. 

The density current a -g^- on the boundary of the domain is a measurable datum. 

We now see how the internal energies can be determined by studying the change 
between a-^- and a-^-. We consider a wavenumber k and two fixed illuminations 

ov av 

tp and tp- Suppose the domain w to be a small ball inside Q of centre z. By 
asymptotic analysis techniques [3 0] , there holds 

(27) [ «(^- d 4)**> 



on 



dv dt 

^ :a(z)Vu£(z) • Vuf(z) + \uj\kc q (z)aq(z)u%(z)ut(z) + o(|w|), 



c a {z)a + d 

where \uj\ denotes the Lebesgue measure of ui and da the integration with respect 
to the surface area. Since the left hand side is known, by choosing different values 
for a the internal power density data 

Ef{z) = a(z)Vut(z) ■ Vufiz), ef(z) = q(z)u*(z)ut(z) 

are recovered for every z £ f2', where fi' <e O is the set of all possible centres where 
the ultrasonic beams are focused. Note that the "polarized" data with ip ^ ip is 
available for a fixed k, but this argument does not allow the reconstruction of the 
cross-frequency data 

Etf(z) = a{z)Vul{z) ■ Vufiz), etf(z) = q(z) u£(*)«f (*). 

We thus chose not to use cross-frequency data for the reconstruction, in contrast 
to[l]. 
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Let us now represent the formulation of Problem l2.4l Given a set of measurement 
K x {ifii} we consider the unique solution u\ £ if 1 (17) to the problem 



(28) 



-div(a Vw|) — kqu l k = 



in £1. 



We define the internal data by 



(29) 



on dQ. 



FP 



For simplicity, we denote eu = {eVr)ij and := and similarly for E. The 
matrices and Ek are to be considered as known matrix-valued functions. We 
study the inverse problem of determining the parameters a and q in f2' from the 
knowledge in fi' of and E^ with a properly chosen set of measurements. Note 
that this reconstruction problem is slightly different to the one studied in [5] , where 
the full matrices 

^kik± ' ' ' &kikM E ki_ki ' ' ' 

E = 



are supposed to be known, with K = {k\, . . 
only the diagonal blocks are measurable. 



&m}- In our case, we suppose that 



5.2. Reconstruction Algorithm. This subsection is devoted to the proofs of the 
exact formulae given in Subsection 12.21 

The new 



5.2.1. Reconstruction Formula for a/q. We restate here Theorem 
contribution of this work is the proof of the bounds ^ . 

Theorem 12.61 Take p, r > and K x {tpi} £ V(Cl';p,r). Suppose that 

© \Ef{x)\,\4{x)\<b, x£Q', 

for some b > 0. Take x £ Q' and k as in Definition 
C = C(A,/?i,b) > such that 

tr(e k )tr(E k )-tr(e k E k ) 



Then there exists 



2„4 



-0) > Cp z r 



and a/q is given in terms of the data by 

/TTnl ivjf U ( \M 2 a o tr ( e fc) tr ( E k) ~ tT ( e k E k) ■ 
(US) |V(efe/tr(efe))| 2 - = 2 in x. 

1 ZT \ e k) 

Proof. Condition ()PSM l|l implies tr(e^) > in Q', and so we may divide by tv(e k ). 
Following the argument given in the proof of Proposition 3.3 in (4j we obtain 

ly7f /. i nn,2 a ti(e k )tv(E- k ) -tr{e k E k ) 
|V(e s /tr<e s ))| a - = 2 ^ . 

We now prove ((9|). For cleanliness of notation we shall denote several positive 
constants depending onA and b simply by c = c(A,6) > 0. Until the end of the 
proof, all the functions will be considered evaluated in x. We equip the space 
of real symmetric matrices with the Hilbert-Schmidt scalar product defined by 
(A, B) = tr(AB). We claim that 



(30) 



2^3, 3 



1 — cos ( 
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As a consequence of this inequality, which we shall prove later, we get 
(31) Jl-cmO e - k , E - k \\E- k \\>c- 



1 



since by flPSM 2ft there holds 

, 2 



^2,2^3,3 



(V«l • Vn? 



> cr z 



Note that, from the definition of and E^, we have 



ef ePJ and 

fc fc 



< 



E^E^ 3 . In particular, 



(32) 



i/E( e ?) a = te (e S ), 11^11 = JEK)^^ 

V *jj V *j 



Combining ([3~T|) and (|32|) we obtain 



tr(ea)tr(g £ )-tr(e^) > ta^)(l - cos0 es , Bj[ ) 



tr(e,) 5 



> 



tr(e,) 2 
:AVtr(e s ) 



" tr( e£ )2 t r(^)(l + lg}) 2 ' 
which gives the desired result since the denominator is bounded by a constant 
depending on b, and tr(e^) > flip 2 ■ 

Let us now turn to the proof of ([50)) . We use the notation g = Ej./ \\Ej.\\ — 
e- k j ||e s ||. By © and ([32]) we have 

2 



J^Z' E- 



< 



(4 3 )' 



+ 



\E 



fell 2,2 £,3,3 



2.3 

e fe 



E- E- 



E- h 



l_ 2,2 3,3 
2 e fc e fe 



2.2^3,3 
" e fe E k 



4" 



\E 



k\\ 2,3 
fe 



I Si 



152,3 1 



2,2 I I 
-efe 1^3,3 1 



EdE?> 3 \g 2 , 2 \ 



<c||S s ||(l- 



tr(^fe) 



COS V e - Er, 



since 1 52,3 1 



153,3 1 



1^2,2 1 < C|jg|| < C^/l - COs0 e -, Eji - 



□ 



5.2.2. Reconstruction Formula for q. In this subsection we derive the reconstruction 
formula for a given in Theorem 12.71 It is based on the knowledge of the ratio 
G := a/q, that can be computed via the formula (j 10|) . Since this reconstruction 
involves the derivative of the data, we need a stable way to obtain q from G. 

The following lemma reviews the derivatives and the trace of products of func- 
tions in Sobolev Spaces. 

Lemma 5.1. Take u,v € if^fi; C) D L°°(0; C). Then uv G H X (Q,;€), V(m/) = 
uVd + wVm and 7(111;) = 7(11)7(1;). // in addition v > c > almost everywhere 
then u/v £ i? 1 (J7;C), S7(u/v) = Vu/w — uVs/w 2 and j(u/v) — r y(u)/ r y[v). 

Proof. It easily follows from [HI Section 5.5, Theorem 1] and [19l Lemma 7.5]. □ 
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We need the following preliminary result, which follows from a straightforward 
calculation by Lemma |5~T1 and Proposition 13.31 The proof is left to the reader. 

Then u l k u j k G H\n) is the 
in SI, 



We are now in a position to prove the reconstruction formula for q given in 
Theorem 12 .71 Recall that if we consider a proper set of measurements K x {pi} we 
have that 

(34) tr(e) = ^ e l k { > c> in 0', 

for some c > 0. 

Theorem [2771 Let Kx{ipi} be a proper set of measurements. Suppose q G H 1 ^). 
Then logq is the unique solution to the problem 

f -div(Gtr(e) Vu) = -div (G V (tr(e))) + 2J2 k , t {E% - kef) in SI', 
\ u — logq\d(i> on dCl'. 

Proof. By Proposition |33] we infer that u k ,u{ G if 1 (SI) n L°°(f2). Therefore, by 
Lemma Owe get that u\u j k G H 1 ^) n L°°(Q). Since q € H\n) n L°°(fi) we 
obtain that e k e 1 (SI) . Hence by Lemma 7.5 in [TH] we have 

VK<) = V (e«/«) = (Vejf - e%VQogq))/q, 

with log g € if 1 (SI). As a result, in view of Lemma I5~2l it is immediate to show that 
for every v G I? (SI) 

{2ke l l - 2EfA (v) = -div(G Ve^')(«) + div fce^ V(logg)) (w), 

whence the result follows by summing all these equations with i = j. □ 

5.3. Numerical Experiments. In this section we shall do some numerical simu- 
lations of the reconstruction algorithm discussed in the last section. FreeFem++ 
has been used. The exact formulae given in Theorem 12.61 and Theorem 12.71 will be 
used to image the electromagnetic parameters a and q. 

In both cases, the construction of proper sets of measurements by means of the 
multi-frequency approach turns out to be effective. Further, in two dimensions, the 
reconstruction procedure gives better results than the one described in [3]: with 
about one seventh of the available data, the reconstruction errors are about half of 
the previous ones. 

5.3.1. Two- Dimensional Example. Since the two-dimensional case has been tested 
thoroughly in [4], we have decided to study the same example in order to be able 
to make a comparison. There are two main differences. First, the formula for the 
reconstruction of q is not the same. Second, in our case the available data is smaller 
since we do not use the cross-frequency data. 



Lemma 5.2. Let K x {ifi} be a set of measurements, 
unique solution to the problem 



(33) 

for every k,i,j. 



-div (oV(«))=2fte^-2E^ 
u = ipiifj on <9S7, 
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Figure 2. The reference parameters. 



2.3 




(a) The coefficient a. 



(b) The coefficient q. 



FIGURE 3. The reconstructed parameters. 




(a) The coefficient a* . 



(b) The coefficient q* . 



Let fl — B(Q, 1) be the unit disk. We use a uniform mesh of the disk with about 
3000 triangles and 1600 vertices. The coefficients are given by 







in B, 


'2 


in B, 


■-I 




in C, 
in E, 


1.8 
1.2 


in C, 
in E, 






otherwise, 


/I 


otherwise 



The set B is the rectangle with diagonal (0, 0.4) — (0.3, 0.5). The set C is the area 
delimited by the curve t i-> (0.3 + p(t) cos(i), — 0.2 + p{t) sin(t)), where p{t) = 
(20 + 3sin(5t) - 2sin(15t) + sin(25<))/100. The set E is the ellipse of centre 
(—0.3, 0.1), with vertical axis of length 0.3 and horizontal axis of length 0.2. These 
parameters represent three inclusions of different contrast in a homogeneous back- 
ground medium. This is the typical practical situation, since cancerous tissues have 
typically higher values in the parameters [28] . The coefficients a and q are shown 
in Figure [2l 

Let f2' = B(0, 0.8) be the subdomain where the internal energies are constructed 
(see Section ISTTj) . We also suppose that a and q are known in ft \ ft'. In view of 
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Table 2. A comparison between the numerical experiments car- 
ried out in [4] and in this work. 





Numerical experiments 

in a 


Numerical experiments 
in this work 


Illuminations 


l,:El, £2 


xi + 2,2:2 + 2 


Frequencies K 


1,3,7 


1,3,7 


Number of energies 


81 


12 


||o-o*|| 2 


3.5 • lCT 1 


1.6 • ICT 1 


Il9-<2*|| 2 


1.5 • ict 1 


0.8 • 10- 1 



Theorem l4.131 we choose the illuminations xi+2 and X2+2 and K = {1, 3, 7}. In [4], 
the boundary conditions 1, xi, x% and the same K were chosen: these illuminations 
satisfy the hypotheses of Theorem 14.131 

Let a* and q* denote the approximated coefficients. We first reconstruct G = 
a* /q* in 1 G O' by means of the formula [TU] 

a* , , = 2 tr(e fc ) tr(E k ) - tr(e k E k ) 

1* \V(e k /tT(e k ))\ 2 2 tv(e k ) 2 W ' 

averaging over all the k € K such that the denominator of the right hand side is 
bigger than 10 -2 . Since for every x S ft' the set of such k is not empty, the chosen 
set of measurements turns out to be proper in ft'. 
Then, we use Theorem 12. 71 to image logg* e iJ^ft'): 

2 

-div(Gtr(e) V(log<f)) - -div(GV (tr(e))) + 2 £ £ - fee") in O'. 

i-eK i=i 

Finally, a* is given by a* = Gg* , which, in absence of noise, gives a good approxi- 
mation. The reconstructed coefficients are shown in Figure [31 

In Table [5] we compare these findings with the numerical experiments performed 
in [1] . Even if the non-availability of the cross- frequency data makes the number 
of measurable internal energies much smaller (about one seventh) , the L 2 norms of 
the errors a — a* and q — q* in this work are about half the corresponding norms 
obtained in [3]. 

5.3.2. Three- Dimensional Example. Take ft = B(0, 1) and ft' = 5(0,0.8). We 
use a mesh with about 30000 tetrahedra and 6000 vertices. For simplicity, we 
choose a = 1 and q = 1 + 0.8 Xb(o,o.3) an d are interested in imaging the elec- 
tric permittivity q. In view of Theorem 14.141 we choose the set of measurements 
{1, 3, 7} x {x\ + 2, X2 + 2}, which a posteriori turns out to be proper in ft'. The same 
reconstruction procedure described in the two-dimensional example is used. The 
reconstruction error is \\q — q*\\ 2 ~ 1.3 • 10 _1 . A comparison between the reference 
and the reconstructed parameters is shown in Figure 01 
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Figure 4. A section of the electric permittivity q in Q'. 




(a) The reference parameter. (b) The reconstructed 

parameter. 
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